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Abstract 

The asymptotic behaviour of vacuum Bianchi models of class A near the initial sin- 
gularity is studied, in an effort to confirm the standard picture arising from heuristic and 
numerical approaches by mathematical proofs. It is shown that for solutions of types 
other than VIII and IX the singularity is velocity dominated and that the Kretschmann 
scalar is unbounded there, except in the explicitly known cases where the spacetime can 
be smoothly extended through a Cauchy horizon. For types VIII and IX it is shown that 
there are at most two possibilities for the evolution. When the first possibility is realized, 
and if the spacetime is not one of the explicitly known solutions which can be smoothly 
extended through a Cauchy horizon, then there are infinitely many oscillations near the 
singularity and the Kretschmann scalar is unbounded there. The second possibility re- 
mains mysterious and it is left open whether it ever occurs. It is also shown that any 
finite sequence of distinct points generated by iterating the Belinskii-Khalatnikov-Lifschitz 
mapping can be realized approximately by a solution of the vacuum Einstein equations of 
Bianchi type IX. 



1. Introduction 

Solutions of the vacuum Einstein equations with an SU(2) isometry group acting on 
spacelike hypersurfaces have been studied extensively over the past twenty-five years or 
more. This class of spacetimes was christened the 'Mixmaster model' by Misner. A useful 
point of entry into the literature on the subject is the book of Hobill et. al. [11], where 
many of the contributions are devoted to the Mixmaster model. Almost all the work which 
has been done in this area is either heuristic in nature or based on numerical calculations. 
There are very few rigorous results. Although the numerical and heuristic approaches have 
led to differences of opinion, it seems that by now a consensus has developed concerning 
various features of the evolution of these spacetimes. This will be referred to in the following 
as the 'standard picture'. This paper is an investigation of what aspects of this standard 
picture can be supported by rigorous theorems. 

There are different motivations which explain the amount of effort which has been 
put into understanding this very special class of spacetimes. One of these is the desire 
to understand the nature of spacetime singularities and in particular to find out whether 
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the geometry near a singularity admits a simple description or whether it is intrinsically 
complicated. The interest of the Mixmaster model is that it is a case where although the 
setting is relatively simple (the Einstein equations reduce to ordinary differential equa- 
tions) the behaviour of the solutions seems to be very complicated. Another motivation 
for studying the Mixmaster model is provided by the idea of Belinskii, Khalatnikov and 
Lifschitz [3] that it could provide an approximate description for very general spacetime 
singularities. This suggestion that very general spacetime singularities can be described 
approximately by a relatively simple model, if true, is only useful insofar as the model 
itself is understood. 

In this paper the ordinary differential equations which describe the Mixmaster model 
will be looked at from the point of view of the theory of dynamical systems. This approach 
has been discussed extensively in the book of Bogoyavlensky [5] . This book contains many 
interesting ideas, but does not contain any general theorems on the Mixmaster model, the 
end result being an impressionistic informal description as to what the behaviour of generic 
solutions should be. Another way of writing the equations for the Mixmaster model as a 
dynamical system has been presented by Wainwright and Hsu [22] . Their equations cover 
a wider class of spacetimes which are spatially homogeneous. In particular they include 
vacuum spacetimes of Bianchi types of class A, i.e. Bianchi types I, II, VIo, VIIo, VIII and 
IX. This means considering spacetimes with three-dimensional isometry groups more gen- 
eral than just SU (2). The type IX spacetimes constitute the Mixmaster model. According 
to the standard picture general spacetimes of type VIII have a complicated singularity 
similar to that of general spacetimes of type IX while the singularities in solutions of the 
other types are much simpler. It turns out to be very helpful to study all class A models 
in a unified way rather that trying to handle type IX in isolation. In [22] various general 
properties of this system were determined, and in particular the nature of its critical points 
was investigated in detail. The work of Wainwright and Hsu is the starting point for the 
investigations of the present paper. 

The results of the paper can be summarized as follows. Any solution of the vacuum 
Einstein equations with a Bianchi symmetry of class A can be classified, according to its 
behaviour near the initial singularity, into one of three types. These will be referred to 
as 'standard convergent', 'standard oscillatory' and 'anomalous'. (For the precise meaning 
of this terminology, see the last paragraph of Section 4). Solutions of types I, II, VIo 
and VIIo belong to the standard convergent type. This is also true of the NUT solutions 
of types VIII and IX, discussed further below. All other solutions of types VIII and IX 
belong to either the standard oscillatory or the anomalous type. If it could be shown that 
the anomalous type never occurs then this would provide a strong confirmation of what 
was referred to above as the standard picture. However, the results of this paper do not 
suffice to obtain that conclusion; in fact they do not even suffice to show that the standard 
oscillatory type ever occurs. 

The detailed conclusions which can be made about the spacetimes of the different types 
just introduced concern on the one hand whether the singularity is a curvature singularity 
and on the other hand precise statements about the convergent or oscillatory nature of 
the spacetimes near the initial singularity. Those vacuum Bianchi models of class A which 
admit an extension through a smooth Cauchy horizon have been determined explicitly in 
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[6] . Here it is shown that in any vacuum spacetime of Bianchi class A which does not admit 
such an extension, and which is not anomalous, the Kretschmann scalar R a p 1 sR a P lS is 
unbounded in a neighbourhood of the initial singularity. This can be seen as a result 
on the non-existence of 'intermediate singularities' in this class of spacetimes. In general 
this result is weakened by the fact that we have no control over how many anomalous 
spacetimes exist. On the other hand, for types other than VIII and IX it provides optimal 
information. As for the question of convergent or oscillatory behaviour, it is shown that 
in the standard convergent case the singularity is velocity dominated. On the other hand, 
in the standard oscillatory case it is shown that naturally defined geometrical quantities 
undergo infinitely many oscillations as the singularity is approached. 

One aspect of the standard picture is that general solutions of the Mixmaster model 
are supposed to be approximated in some sense near the singularity by a mapping of 
the circle to itself, the BKL (Belinskii-Khalatnikov-Lifschitz) mapping. The results on 
oscillations mentioned up to now say nothing about this. However it can be proved to be 
true in a weak sense. The statement is that, given any finite sequence generated by the 
BKL map, there exists a solution of the Einstein vacuum equations of type IX (or indeed 
of type VIII) which reproduces this sequence with any desired degree of accuracy. 

The paper is organized as follows. Section 2 presents some parts of the theory of 
dynamical systems which will be necessary in the analysis. In Section 3 the system of 
Wainwright and Hsu is recalled. Various results on the asymptotic behaviour of a solution 
of the Wainwright-Hsu system as the singularity is approached are proved in Section 4. 
In Section 5 these results are used to obtain conclusions about curvature singularities 
and the convergent or oscillatory behaviour of the corresponding spacetime. Moreover, it 
contains the precise formulations of the results which have been presented informally in 
the preceding paragraphs. 



2. Background on dynamical systems 

The purpose of this section is to collect together some facts on dynamical systems 
which are useful in the analysis of the dynamics of Bianchi models. For a more general 
introduction to relevant aspects of the theory the reader is referred to [21]. Consider the 
system of ordinary differential equations: 

d i = m (2.D 

where / : R n — > R n is a C°° mapping. (The choice of the differentiability class C°° is not 
essential in what follows; it suffices for the applications in this paper.) The corresponding 
local flow is the mapping F defined by the condition that F(t, x) is the value at time t of 
the solution of (2.1) which takes on the value x at time t = 0, provided the solution exists 
that long. A critical point of (2.1) is a point x with f(x) = 0. An important element in 
the analysis of the global properties of solutions of (2.1) is the study of the local behaviour 
of solutions near critical points. This can be described using the concept of topological 
equivalence. 

Definition If xq is a critical point of (2.1) and yo is a critical point of the system dy/dt = 
g(y) with local flow G(t,y), then the two systems are said to be topologically equivalent 
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near the points xq and yo respectively if there exists a homeomorphism <fi of an open 
neighbourhood U of xq onto an open neighbourhood V of yo, with 4>(xq) = yo such that 
whenever F(t,x) is defined for some x G U, G(t,<f)(x)) is also defined, and is equal to 
(f>(F(t,x)). 

In other words, when expressed in appropriate local coordinates the solutions of the two 
systems look identical. A critical point xq of the system (2.1) is called hyperbolic if the 
derivative Df(xo) of / at that point (which defines the linearization of the system about 
that point) has no eigenvalues which are purely imaginary. A fundamental result on 
hyperbolic critical points is the Hartman-Grobman theorem ([9], p. 244): 

Theorem 2.1 If xo is a hyperbolic critical point of the system (2.1) then the system is 
topologically equivalent near xq to the linearized system dx/dt = Df(xo)x near the origin. 

Note that the homeomorphism whose existence is required in the definition of topological 
equivalence cannot in general be chosen to be a diffeomorphism. However this fact will 
play no role in what follows. Linear systems which have no purely imaginary eigenvalues 
can be classified up to topological equivalence ([2], p. 48). Let n + and n_ be the number 
of eigenvalues with positive and negative real parts respectively. Two linear systems of 
ODE without purely imaginary eigenvalues are topologically equivalent near the origin if 
and only if the corresponding values of n + and n_ are equal. (In fact in this case the 
equivalence is global, i.e. U and V can be taken to be all of R n .) Thus up to topological 
equivalence, the only linear systems on R n with no purely imaginary eigenvalues are given 
by the systems on R n + x R n - defined by dy/dt = y, dz/dt = —z, with y G R n+ and 
z G R n ~ . A system of ODE of this type is known as a standard saddle and there are only 
Ti+1 possibilities in R n . If n + or n_ vanishes then there is not a saddle in the usual 
sense, but rather a source or sink. Nevertheless, we include this case in the definition of 
a standard saddle. Combining this discussion with the Hartman-Grobman theorem shows 
that in a neighbourhood of any hyperbolic critical point the system (2.1) is topologically 
equivalent to a standard saddle. 

What can be said in the case of a critical point xq which is not hyperbolic? Let E+, 
Eq and E_ be the spaces spanned by those generalized eigenvectors of Df(xo) whose real 
parts are positive, zero and negative respectively. These are called the unstable, centre and 
stable subspaces respectively. In general R n is the direct sum of these three subspaces and 
the hyperbolic case is that where the centre subspace reduces to zero. Associated to these 
subspaces are locally invariant manifolds. A submanifold M of an open neighbourhood U 
of x is called locally invariant if whenever x G M and F(t, x) G U then F(t, x) G M. 

Definition A (local) stable, centre or unstable manifold of the critical point x is a C 1 
submanifold of an open neighbourhood U of xq which is locally invariant, which contains 
xq, and whose tangent space at xq is the stable, centre or unstable subspace respectively. 

From the topological classification of hyperbolic fixed points discussed above it is clear 
that a hyperbolic fixed point has stable and unstable manifolds which are locally unique, 
i.e. the intersection of any such manifold with a sufficiently small neighbourhood of the 
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critical point is unique. A generalization of this to non-hyperbolic critical points is given 
by the centre manifold theorem [1]: 

Theorem 2.2 Let xq be a critical point of (2.1) and suppose that / is C°° . Then there 
exist stable and unstable manifolds of class C°° and a centre manifold of class C k for any 
finite k. The stable and unstable manifolds are locally unique. 

The centre manifold need not be C°° in general; the C k centre manifolds may shrink as k 
increases. The centre manifold need also not be locally unique. (See [23], p. 210, for an 
example). Despite the non- uniqueness of the centre manifold, it can be used to formulate 
a generalization of the Hartman-Grobman theorem to non-hyperbolic critical points. This 
is the reduction theorem of Shoshitaishvili ([18], [19]). For a detailed proof in English see 
[13]. 

Theorem 2.3 Let x be a critical point of (2.1). Then the system is topologically equiva- 
lent near x$ to the Cartesian product of a standard saddle with the restriction of the flow 
to any centre manifold. 

Suppose for a moment that / is such that all solutions of (2.1) can be extended so as to 
be defined for all real values of t. Then an ct-limit point of a solution x(t) of the system 
(2.1) is a point x* such that there exists a sequence t n with t n — > — oo and x(t n ) — > x*. 
The a-limit set of the solution is the set of all a-limit points. The concepts of w-limit 
point and cu-limit set are defined analogously by replacing — oo by oo. If not all solutions 
of (2.1) extend to global in time solutions it is possible to rescale / with a positive function 
so that all solutions of the rescaled system do exist globally in time. This is discussed by 
Wainwright and Hsu[22], section 3. The notions of a- and cu-limit points can be applied to 
the rescaled system. The image of a solution of the original system in R n is also the image 
of a solution of the rescaled system. It is natural to define the a- and ui- limit points of a 
solution x(t) of the original system to be those of the corresponding solution of the rescaled 
system. This definition can be rephrased in terms of the original system. A solution x(t) 
of (2.1) defined on some time interval is said to be maximally extended if it cannot be 
extended to a solution on any strictly longer interval. If a maximally extended solution 
xit) is defined on the interval (£_,£+), then the ct-limit and cu-limit points of the solution 
are limits of sequences of the form x(t n ) where t n — > t- or t n — > £+. This is consistent 
with the previous definition in the case that the solution is defined for all real values of t. 

Some standard properties of a-limit sets will now be listed (cf. [23], p. 46) Corre- 
sponding properties hold for u;-limit sets. The a-limit set of any solution is closed. It 
consists of a union of images of solutions of the ODE. It is clear that it must also contain 
the ct-limit sets of these solutions. If the solution stays in a compact set for all t < t 
then the a-limit set is connected. Monotone functions are a useful tool for locating a- and 
w-limit sets, as shown by the following simple lemma. 

Lemma 2.1 Let U be an open subset of R n and let F be a continuous function on U such 
that F(f(t)) is strictly monotone for any solution f(t) of (2.1) as long as f(t) is in in U. 
Then no solution of (2.1) whose image is contained in U has an a- or w-limit point in U . 
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Proof Suppose that p G U is an o-limit point of a solution f(t) whose image is contained in 
U. Then F(f(t)), being strictly monotone, must have a limit, possibly infinite, as t — > £_. 
On the other hand, there is a sequence t n with t n — > t- such that F(f(t n )) coverges to 
F(p). Hence F(f(t)) converges to F(p). Thus F is constant on the o-limit set of f(t). 
There exists a solution / (t) which passes through p and is entirely contained in the o-limit 
set of f(t). It follows from the above that F is constant along f(t) : contradicting the 
property of strict monotonicity. The argument for the w-limit set is similar. 

3. The equations of Wainwright and Hsu 

There are many popular ways of writing the equations for Bianchi models. The 
analysis which follows uses a form of the equations for models of class A due to Wainwright 
and Hsu [22]. This is natural since it is in a sense an extension of their approach. The 
results of the present paper concern vacuum spacetimes, whereas the Wainwright-Hsu 
system describes a perfect fluid with linear equation of state p = (7 — l)p. However, 
vacuum models are described by the restriction to the submanifold defined by the vacuum 
Hamiltonian constraint. The latter is independent of 7 and so the choice of 7 is immaterial 
for our purposes. However, the explicit form of the equations off the constraint hypersurface 
is more or less complicated according to the value of 7. In order to take advantage of this, 
7 = 2/3 is chosen in the following. With this choice, the equations are: 

N[ = (q-4E+)N l 

N' 2 = ( 9 + 2£ + + 2 v / 3£_)iV 2 

= (q + 2E+ -2V3E_)AT 3 (3.1) 
£' + = -(2-g)E + -3S+ 
£'_ = -(2 - q)E_ - 3S- 

where 

q = 2(E 2 + + J: 2 _) 

S+ = |[(iV 2 - N 3 ) 2 - ^(2^ -N 2 - N 3 )] (3.2) 
S_ = ^(N 3 - iV 2 )(iV 1 -N 2 - N 3 ) 

and a prime denotes a derivative with respect to a certain time coordinate r. The (vacuum) 
Hamiltonian constraint is: 

£+ + T?_ + § [N 2 + N 2 + N 2 - 2(N±N 2 + N 2 N 3 + N3N1)] = 1 (3.3) 

A solution of the vacuum Einstein equations with a Bianchi symmetry of class A is de- 
scribed by a solution of (3.1) with initial data satisfying (3.3). Then of course the whole 
solution lies in the submanifold defined by (3.3). If t is a Gaussian time coordinate based 
on one of the orbits of the group action defining the symmetry, and if trk(t), denotes the 
mean curvature of these orbits, then the time coordinate r in (3.1) is defined by the relation 
r(£) = — I f* trk(t')dt' . Since the determinant of the induced metric of the orbits satisfies 
d /dt (log detg) = — 2tr/c, r 3 is proportional to the volume form of the orbits. In other 
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words r represents a length scale related to that volume form. For the precise definition of 
the variables (N\, N 2 , N3, E + , E_) the reader is referred to [22]. The variables N±, N 2 and 
N3 describe the curvature of the spatial slices and, at the same, the Bianchi type. In par- 
ticular they allow all Bianchi types of class A to be included in a single dynamical system. 
When a fluid is present, E + and E_ are related to the shear. In the vacuum case they 
represent the trace-free part of the second fundamental form. All of these variables are 
dimensionless, in the sense that if the spacetime metric is multiplied by a constant they do 
not change. This means that they are candidates for quantities which remain finite when 
a spacetime singularity is approached. The system (2.3) has a three-fold symmetry, which 
is not obvious when the equations are written in this form. It is described in detail in [22] . 
In particular, it cyclically permutes Ni, N 2 and N 3 . 

In discussing the global structure of the spacetimes considered in this paper, we restrict 
to the maximal Cauchy development of data on a spacelike orbit, which is in any case the 
only part of the spacetime described directly by the variables in (3.1). It is well-known 
that a vacuum Bianchi model of class A, and not of type IX, is geodesically complete in one 
time direction which, reversing the time orientation if necessary, we can assume to be the 
future. In the other time direction there is a finite upper bound to the length of all causal 
geodesies. In the case of Bianchi type IX, a bound of this kind holds in both directions as a 
consequence of a result of Lin and Wald[14,15]. In the latter case the variables in (3.1) are 
not defined at the moment of maximum expansion. For the definition of these variables 
involves dividing by the mean curvature tr/c, which vanishes at that moment. Thus a 
maximally extended solution of (3.1) of Bianchi type IX represents only half of a maximal 
Cauchy development, where the model is either always expanding or always contracting. 
Reversing the time orientation if necessary, it can be assumed that it is expanding, so that 
the singularity lies in the past. With these conventions, if a maximal solution of (3.1) is 
defined on the interval (t_,t+) then the behaviour of the corresponding spacetime near 
the singularity is encoded in the behaviour of the solution as r — > r_. For Bianchi types 
other than IX the behaviour as r — > r + encodes the behaviour of the spacetime in the time 
direction where it is complete. For Bianchi type IX it encodes the behaviour near the time 
of maximal expansion. 

The Bianchi type to which a solution of (3.1) corresponds depends on the values of 
N1-N3. If all three are zero the Bianchi type is I. If precisely one is non-zero then it is II. 
If precisely two are non-zero it is either VIo (signs opposite) or VIIo (signs equal). If all 
three are non-zero it is either IX (all signs equal) or VIII (one sign different from the other 
two). The set of points corresponding to any one of the Bianchi types is invariant under 
the flow of the differential equation. A point of this set will be referred to as a point of 
the given Bianchi type. 

There are various monotone functions which can be defined. The function NiN 2 N 3 is 
strictly monotone along any solution of type IX or VIII. (Cf. the function Ai of [22].) For 
(N 1 N 2 N 3 y = 3q(N 1 N 2 N 3 ). This gives the result for q ^ 0. If q = then the Hamiltonian 
constraint shows that either E' + or E'_ is non-zero and this proves the desired result. The 
functions 

Z e = [f 3£ + (at 2 + e Ar 3 )2]/(_ ei v 2 iV3) 
are non-negative functions on the sets of points of type VIq and VIIo with Ni = for 
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e = 1 and e = — 1 respectively. They are strictly decreasing except when E_ = 0. These 
statements are proved in [22]. The function (1 + E + ) 2 is a non-increasing function along 
solutions on the union of points of types I, II, VIo and VIIo with N± = 0. For 

d/dr{{\ + E+) 2 ] = -4(1 + E+) 2 (l -Ti 2 + - E -) ( 3 - 4 ) 

and the vacuum Hamiltonian constraint implies that E 2 ^ + T?_ < 1 for the given Bianchi 
types. This is analogous to a monotone function for solutions of class B given by Hewitt 
and Wainwright [10]. 

In order to study the question of which of the spacetimes described by the above equa- 
tions have curvature singularities, it is useful to express curvature invariants of the space- 
times in terms of the Wainwright-Hsu variables. Consider the example of the Kretschmann 
scalar R^p^sR 01 ^ 5 . Define a dimensionless version of this quantity by 

k = R^sR^/itrk) 4 (3.5) 

Then k can be expressed as a polynomial in the variables AT 1; jV 2 , N 3 , E + and E_. This 
can conveniently be done using the formulae for curvature given by Ellis and MacCallum 
[8]. Here only the expression for solutions of type I will be written out. In that case: 

«=(-§ + |S+ + f Tlf + 2(1*1 + §E 2 _ - ±E + - |) 2 + *E 2 _ (2E + - I) 2 (3.6) 

The fact that it is a positive definite expression is explained by the fact that for a vacuum 
spacetime of Bianchi type I the magnetic part of the Weyl tensor vanishes. 

4. Possible limit sets 

It will now be shown among other things that the ct-limit set of any solution of (3.1)- 
(3.3) of type I, II, VIo or VIIo is a single point of type I or type VIIo- Consider first 
a solution of type I. It corresponds to a critical point of (3.1) and so is its own a-limit 
set, as well as its own cu-limit set. These points form a circle, sometimes known as the 
Kasner ring. Next consider a solution of type II. Using the threefold symmetry of the 
equations, it may be assumed without loss of generality that it satisfies N 2 = N% = 0. 
Then an elementary calculation shows that either E_ = or the ratio (E + — 2)/E_ is 
time independent. Taking account of the Hamiltonian constraint, this reduces the motion 
to motion along a curve. The direction of motion along this curve is controlled by the 
monotone function (1 + E + ) 2 . It follows that the a-limit set consists of one point on the 
Kasner ring and the w-limit set of another. There are three points on the Kasner ring 
which play a special role. They are denoted by T\ , T 2 and T3 in [22] and have coordinates 
(— 1, 0), (1/2, ±-\/3/2). They divide the Kasner ring into three equal parts. The a-limit set 
of a solution of type II with N\ ^ lies on the longer of the two open arcs with endpoints 
T 2 and T3, while the w-limit set lies on the shorter of these arcs. The points T\, T 2 and T3 
are permuted cyclically by the threefold symmetry, which shows what happens for type II 
solutions with N 2 7^ or ./V3 7^ 0. Associating the ct-limit set of a type II solution with 
the o;-limit set of that solution defines a mapping from the Kasner ring with the three 
exceptional points removed to the Kasner ring. This is a realization of the famous BKL 
(Belinskii-Khalatnikov-Lifshitz) mapping. 
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Consider next solutions of type VI or VII . Using the threefold symmetry, it can 
be assumed without loss of generality that N\ = 0. Equation (3.4) and Lemma 2.1 imply 
that for any a-limit point £+ + = 1 or E + = — 1. For these solutions the Hamiltonian 
constraint reduces to: 

Z\ + Z 2 _ + l(N 2 - N 3 ) 2 = 1 (4.1) 

and so if E + = — 1 then + T?_ = 1 anyway. Moreover N 2 = N 3 . This means in 
particular that there are no a-limit points of type VIo- The a-limit set must contain a 
solution of the equations, which satisfies these conditions everywhere. The second and 
third equations of (3.1) then give either N 2 = N3 = or £_ = 0, whence = ±1. 
Where E + = — 1, E + + 1 = and since this function is no n- negative and non-increasing 
along solutions, it follows that any solution which has a point of this type as an a-limit 
point must be a time independent solution. Apart from these time-independent solutions, 
the remaining possibilities are that either the a-limit set contains a point with E+ = 1, 
E_ = and N 2 = N3, or that the entire a-limit set is contained in the Kasner ring. The 
points in the first of these two cases constitute the image of a whole solution and so if 
one of them belongs to the a-limit set of a solution, they all do. In particular, this would 
mean that the solution would not remain in a compact set. It will now be seen that this is 
impossible. In the Bianchi VIo case this is simple, since the Hamiltonian constraint shows 
that the set of Bianchi VIo points is compact. The Bianchi VIIo case requires a little more 
work. Consider the behaviour of the non-increasing function Z_\ along any solution of 
type VIIo- Either it is identically zero or or it is bounded below by a positive constant 
for all t < to and some fixed tq. In the first case N 2 = N3, E_ = and E + = ±1. Then 
N2 = 2(1 + T. + )N 2 and the solution remains in a compact set as r — > — r_. In the second 
case, for any fixed To, 

N 2 N 3 < C[p?_ + (N 2 - N 3 f] (4.2) 

holds for all r < tq. Combining this with the Hamiltonian constraint (4.1) proves that 
the solution remains in a compact set as r — > — r_. The fact that solutions of these types 
remain in a compact set for all sufficiently negative times also proves that their a-limit 
sets are non-empty. 

The points Ti, T 2 and T3 of the Kasner ring correspond to spacetimes obtained by 
making identifications in a subset of Minkowski space. In particular they are flat. In fact 
they are the only flat Kasner solutions, since from (3.6) the Kretschmann scalar is non- 
zero at all other points of the Kasner ring. The points of type of type VIIo which satisfy 
E_|_ = —1 and E_ = represent the same spacetime as the point T\. The reason for this 
is that the Wainwright-Hsu variables depend not only on the spacetime but on a choice of 
frame. The flat Kasner solutions admit a frame of type I and a one-parameter family of 
frames of type VIIo, which can all be used. Thus one spacetime is represented by different 
points in R 5 . These type VIIo solutions and those obtained from them by applying the 
threefold symmetry will be referred to as the flat type VIIo solutions. This terminology 
will be justified later, when it is shown that they are the only flat solutions of type VIIo- 

It has now been shown that, except for the flat solutions of type VIIo, the a-limit 
set of any solution of type I, II, VIo or VIIo is contained in the Kasner ring. It will now 
be shown that it consists of a single point of the Kasner ring. To see this, consider the 
behaviour of the function (1 + E + ) 2 along the given solution. It cannot take two different 
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values on the a-limit set, since then it could not be monotonic along the solution. Thus 
the a-limit set of any given solution consists of points on the Kasner ring where E + takes 
on the same value. But for any given value of E + there is only one point, in which case the 
desired result follows, or there are two. In the latter case it follows from the connectedness 
of the a-limit set that only one of the two can be contained in it. The following theorem 
has now been proved: 

Theorem 4.1 The a-limit set of a solution of (3.1)-(3.3) of type I, II, VIo or VIIo is a 
single point of type I or a flat point of type VIIo, the latter only being possible if the 
solution is time independent. 

Note that all the arguments used to prove this were elementary, and that the reduction 
theorem (Theorem 2.3) was not used. In studying the a-limit sets of solutions of types VIII 
and IX the heavier machinery is required. Note first that, due to the strictly monotone 
function for these Bianchi types presented in Section 3, the a-limit set of any solution of 
type VIII or IX is contained in the set NiN 2 Ns = 0, which consists of points of the other, 
simpler, Bianchi types. It follows that the a-limit set, if non-empty, contains the image 
of a solution of one of these types. It then also contains the a-limit set of that solution. 
Thus it follows from Theorem 4.1 that: 

Theorem 4.2 The a-limit set of a solution of (3.1)-(3.3) of type VIII or IX is either empty 
or contains a point of type I or contains a flat point of type VIIo- 

In the rest of this section it will be shown that, except for a small set of well-understood 
solutions, if the a-limit set of a solution of type VIII or IX is non-empty, then it either 
contains a non-empty open subset of flat points of type VIIo or a non-flat point of the 
Kasner ring and that in the latter case it contains more than one non-flat point of the 
Kasner ring. 

Theorem 4.3 The a-limit set of a solution of (3.1)-(3.3) of type VIII or IX has at least 
one of the following properties: 

(i) it is empty 

(ii) it contains at least two distinct points of type I, at least one of which is non-flat 

(iii) it consists of flat points of type I or VIIo and contains a non-empty open subset of the 
set of points of type VIIo- 

(iv) it is a flat point of type I or VIIo and the solution is contained in the unstable manifold 
of that point 

Proof The essential point is to apply the reduction theorem (Theorem 2.3) to the critical 
points of the system. In this context it is important to note that the Kasner ring is a 
centre manifold for each of the non-flat points contained in it, while the manifolds Ni = 0, 
N2 = and N3 = are centre manifolds for the flat solutions of type I and VIIo which 
they contain. This follows from the computation of the eigenvalues and eigenspaces of 
the critical points carried out by Wainwright and Hsu [22] and the obvious fact that the 
manifolds in question are invariant manifolds. 

Consider now a solution of type VIII or IX with non-empty a-limit set. First it will be 
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shown that if its a-limit set contains one non-flat point p of type I, it contains at least two 
such points. The Kasner ring, which is a centre manifold for p, consists of critical points. 
Hence, by the reduction theorem, a neighbourhood of p is foliated by invariant manifolds 
of codimension one, each of which meets the Kasner ring in a single point. In each of 
these submanifolds the flow is topologically equivalent to a standard saddle. It follows 
that if the solution converges to p as r — > r_ , it must lie on the unstable manifold of p and 
that otherwise the ct-limit set must contain the stable manifold. Suppose without loss of 
generality that p belongs to the shorter of the two arcs of the Kasner ring joining T 2 to T3. 
Then the analysis of the linearization in [22] shows that the subspaces N 2 = A3 = and 
Ni = are the stable and unstable manifolds respectively. Since the unstable manifold 
contains no points of types VIII or IX, it follows that the a-limit set of the solution under 
consideration contains the stable manifold, which is the image of a solution of type II. 
Since the a-limit set is closed, it follows that it must also contain the image of p under the 
BKL map. Thus it contains at least two distinct points. 

Next the case will be considered where no a-limit point of the solution is a non-flat 
point of type I. Because of the monotone function NiN 2 N 3: all ct-limit points must satisfy 
Ni = 0, N 2 = or A 3 = 0. A a-limit point of type II or type VIo or a non-flat point of 
type I or VIIo is not possible, as follows from Theorem 4.1. Hence the entire a-limit set 
consists of flat points of type I or VIIo- Note also that no solution on the centre manifold 
can approach one of these points as r — > r_ , unless it is the solution which stays at that 
point all the time. Hence, by the reduction theorem, the solution under consideration can 
only approach a point of that type if it lies on the stable manifold of that point. Using 
the calculations of [22] once more, the stable manifold can be identified. After applying 
the threefold symmetry if necessary, it is given by the equations N 2 = A3, £_ = 0. 
The Bianchi type IX solutions satifying these conditions are the Taub-NUT solutions. 
Similarly, solutions of type VIII (or indeed of type II) satisfying these conditions are the 
NUT solutions of those types described by Siklos[20]. It remains to consider the case of a 
solution whose a-limit points are flat points of type I or VIIo, but which does not converge 
to such a point. By the reduction theorem, a solution of this type which has p as an 
a-limit point must repeatedly leave the neighbourhood of p whose existence is asserted 
by the theorem. In particular, it crosses a sphere of any sufficiently small radius about p 
infinitely many times. Thus some a-limit point lies on any such sphere. It follows that the 
a-limit set contains an open subset of the set of flat points of type VIIo- 

The notions of 'standard convergent', 'standard oscillatory' and 'anomalous' mentioned 
in the introduction will now be defined in terms of the four cases of the theorem. A 
solution of type VIII or IX is standard convergent if it satisfies property (iv) of Theorem 
4.3. A solution of any other Bianchi type of class A is always standard convergent. (This 
terminology is justified by Theorem 4.1.) A solution is standard oscillatory if it satisfies 
property (ii) of Theorem 4.3. A solution which is neither standard convergent nor standard 
oscillatory is called anomalous. 

5. Conclusions 

The results of the previous section on the a-limit sets of solutions of (3.1)-(3.3) will 
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now be interpreted in terms of properties of the corresponding spacetimes. The first issue 
to be considered is that of curvature singularities. Here the theorems of the last section 
give complete information for all types other than VIII and IX. For these two types the 
information obtained is less satisfactory. 

Theorem 5.1 Let (M,g) be a vacuum spacetime with a Bianchi symmetry of class A. 
Suppose that the time orientation has been chosen such that the maximal Cauchy devel- 
opment of data on a homogeneous hypersurface is past incomplete, so that there is a past 
singularity. Then at least one of the following holds: 

(i) the Kretschmann scalar is unbounded in a neighbourhood of the initial singularity 

(ii) the maximal Cauchy development can be extended through a smooth Cauchy horizon, 
and then the spacetime is flat or a NUT spacetime 

(iii) the corresponding solution has property (i) or property (iii) of Theorem 4.3 
Proof Assume that property (iii) of the conclusions of this theorem does not hold. Then 
if the solution of (3.1)-(3.3) corresponding to the given spacetime is of type VIII or IX, 
it has property (ii) or (iv) of the conclusions of Theorem 4.3. The solutions which have 
property (iv) of Theorem 4.3 have already been identified; they are the NUT spacetimes. 
It is known that they admit an extension through a smooth Cauchy horizon (cf. [6]). Thus 
it remains to prove the theorem in the case that a solution of type VIII or IX has property 
(ii) of Theorem 4.3 and in the case that the Bianchi type is neither VIII nor IX. From 
Theorems 4.1 and 4.3 it follows that either the spacetime corresponds to a flat point of 
type I or VIIo or that there is a non-flat point of type I in the a-limit set. In the first 
of these cases, it is known that the spacetime can be extended through a smooth Cauchy 
horizon [6]. In the second case the values of the scalar k introduced in Section 2 have a 
cluster point k$ ^ as r — > r_ , where is the value of k at a non-flat a-limit point of the 
solution on the Kasner ring. Since the Kretschmann scalar is equal to K(tr/c) 4 , it follows 
that lim sup R af3lS R^p^s = oo. This completes the proof. 

According the standard picture of the Mixmaster solutions, case (iii) of Theorem 5.1 should 
not occur. If that were the case, a clean result would be obtained on the non-existence of 
'intermediate singularities'. This would mean that every solution possessed either a smooth 
Cauchy horizon or an unbounded curvature invariant near the singularity. Note that this 
conclusion has been proved for solutions of type other than VIII and IX, as can be seen 
from an examination of the above proof. In fact, slightly more has been proved, since it has 
not only been shown that the Kretschmann scalar is unbounded near the singularity, but 
also that it tends to infinity there at the same rate as rate (tr/c) 4 . Computing curvature 
invariants other than the Kretschmann scalar would not lead to an improved result for 
types VIII and IX using the above techniques, since all curvature invariants vanish at the 
flat points. 

The other important point of interpretation concerns the question of convergent or 
oscillatory behaviour near the singularity, and velocity dominance. The notion of a velocity 
dominated singularity was introduced by Eardley, Liang and Sachs [7] and a related notion 
was used by Isenberg and Moncrief [12] (see also [16], [17] for other applications). The 
general idea is that near the singularity the spacetime should be approximated in some 
appropriate sense by a solution of Bianchi type I. In the present context it might be 
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tempting to say that a solution of (3.1)-(3.3) was velocity dominated if its a-limit set 
consisted of a point on the Kasner ring. However, this is too simple, since the flat points 
of type VIIo represent the same spacetimes as the flat points of Bianchi type I. Thus the 
following definition will be used: 

Definition A singularity in a spacetime with a given 3+1-decomposition is called weakly 
velocity dominated if the generalized Kasner exponents Pi(t, x) converge to a limit for each 
fixed x as the singularity is approached and if their limits satisfy T>ip\ = 1. 

This definition requires some further explanation. Suppose that the hypersurfaces of con- 
stant t have nowhere vanishing mean curvature near the singularity. (If this is not the case, 
the definition is deemed to be violated) Let A^ be the eigenvalues of the second fundamen- 
tal form. Then the mean curvature is trk = E^A^ and the generalized Kasner exponents 
are defined to be pi = Ai/(tr/c). In a Kasner solution they are constants and satisfy the 
equation required in the above definition. For the Bianchi spacetimes considered here, 
the definition will be applied with a 3+1-decomposition defined by Gaussian coordinates 
based on a homogeneous hypersurface. Then the generalized Kasner exponents are linear 
combinations of £_ and E_|_ and the equation of the definition becomes + T?_ = 1. 
It is easy to think up alternative definitions. For instance, the condition T^pf = 1 could 
be replaced by the condition pi/(trk) 2 — > 0, where pi are the eigenvalues of the spatial 
Ricci tensor. The resulting definition is a priori stronger than that given above but in all 
the cases where the above definition is shown to be satisfied in the following the stronger 
definition is also satsisfied. 

Returning to the question of oscillatory behaviour, a function F on an interval (r_, r+) 
will be said to have infinitely many oscillations as r — > r_ if there exist two numbers a, b 
with a < b and a sequence r n with r n — > r_ as n — > oo such that for any positive integer 
k, F(r 2 fc-i) < a and F(r 2 k) > b. 

Theorem 5.2 Let (M, g) be a vacuum spacetime with a Bianchi symmetry of class A. 
Suppose that the time orientation has been chosen such that the maximal Cauchy devel- 
opment of data on a homogeneous hypersurface is past incomplete, so that there is a past 
singularity. Then either: 

(i) the singularity is weakly velocity dominated 

(ii) the spacetime is of type VIII or IX and E + or £_ has infinitely many oscillations as 
the singularity is approached. 

(iii) the spacetime is of type VIII or IX and the corresponding solution of (3.1)-(3.3) fails 
to be contained in a compact set as r — > r_. 

Proof If the spacetime is not of type VIII or IX then the result follows immediately from 
Theorem 4.1. Consider a spacetime of type VIII or IX such the corresponding solution 
of (3.1)-(3.3) is contained in a compact set as r — > r_. Then it must have one of the 
properties (ii), (iii) or (iv) of Theorem 4.3. If it has property (iv) of that theorem the 
singularity is weakly velocity dominated. If it has property (ii) it repeatedly comes close 
to two different points of the Kasner ring as r — > r_ and so E + or £_ has infinitely many 
oscillations as the singularity is approached. If it does not belong to any of these cases 
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of Theorem 4.2 then its a-limit set must be a compact subset of the set of flat points of 
type I and VII . In that case E + — > —1 and £_ — > and the solution is weakly velocity 
dominated. 

This result implies in particular that solutions which are neither of type VIII or IX are 
(weakly) velocity dominated. The standard picture indicates that the only solutions of 
types VIII and IX which are velocity dominated are the NUT solutions and that case (iii) 
of Theorem 5.1 would be superfluous. However, it has not been proved that this is the 
case. 

Although the BKL mapping plays an implicit role in these results, they do not provide 
any direct confirmation of the idea, which is part of the standard picture, that the BKL 
mapping represents an approximation to solutions of the Mixmaster model in some sense. 
The next result is a statement which goes in that direction. Recall that the BKL mapping 
is a mapping from the Kasner ring with the three exceptional points T\ , T 2 and T 3 removed 
to the Kasner ring. Starting with a point of the Kasner ring and applying the BKL mapping 
repeatedly produces a sequence of points. If one of these points is one of Ti, T 2 or T 3 then 
the map can no longer be applied and only a finite sequence can be defined. If, on the 
other hand, this never happens, an infinite sequence is obtained. This may or may not be 
periodic. By a finite non-repeating sequence of BKL iterates we mean a finite initial piece 
of the iteration, which never hits the same point twice, whether or not the sequence as a 
whole is finite or infinite for that starting point. In the case that the whole sequence is 
periodic, we need to cut it off after some number of iterations smaller than the period. To 
measure the quality of the approximation the Euclidean metric in the space R 5 where the 
Wainwright-Hsu system is defined will be used. The distance between points x and y in 
this metric will be denoted by d(x, y). 

Theorem 5.3 Let {xi, . . . , x n } be a finite sequence of BKL iterates and e > 0. Then there 
exists a 5 > and a solution of the vacuum Einstein equations of Bianchi type IX which, 
when written in the Wainwright-Hsu variables as /(r), has the following properties: 

1. 5 < e and the balls of radius 5 about each of the points Xi have disjoint closures. 

2. There exists a finite sequence of times T\ > . . . > r n such that /(r^) is contained in the 
open ball of radius 5 about Xi while for n < r < Tj+i the solution does not come closer 
than 5 to any Xj other than Xi or Xi+i. 

Proof As explained in Section 3 each iteration of a point via the BKL map can be repre- 
sented by a solution of Bianchi type II. For i = 1, 2, . . . , n — 1, let Sj be the image of the 
Bianchi II solution which produces Xi+i from Xi. Choose some rj > smaller than e such 
that the closures of the balls Bi of radius rj about the points Xi are disjoint and such that 
on each of these balls the flow has the local product structure of the reduction theorem. 
Now choose 5q < rj such that any point which is a distance less than 8q from each of two 
distinct Sj must be inside one of the balls Bi. That such a Sq exists follows from the fact 
that the intersections of each Sj with the complement of the union of the Bj are disjoint 
compact sets. Let S be the set of all points which are a distance less than Sq from some Sj. 
It is an open set which contains all Xi. For 1 < % < n — 1 let yi be a point of Sj which is a 
distance less than 8q from Xi. For 2 < i < n let Zi be a point of s$_i which a distance less 
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than 5q from Xi. Let 8\ be such that the ball of radius 8\ about y\ is contained in the ball of 
radius 8q about x\. There exists a time At\ such that the solution with /(r) = z<i satisfies 
/(t + Ari) = yi. Choose 62 > such that if d{x,Z2) < £2 then <i(F(Ari, x), yi) < Si. 
Here F denotes the local flow of the differential equation. By reducing the size of 62 if 
necessary, it can be ensured that if d(x, z 2 ) < £2, then F(t,x) G S^-E^ for all r < Ati, 
where B\ is the complement of Bi UBi+i. Now choose 82 such any x with the correct signs 
of Ni, N 2 and N3 which lies in the open ball of radius 82 about y2 passes through the 
open ball of radius 62 about Z2 and stays in £ as long as it stays in the ball of radius 8q 
about X2- The existence of a number 82 with this property follows from the local saddle 
point structure of the flow near the point X2- In a similar way we can recursively define 
ATj for i = 2, . . . , n — 1, ti for i = 3, . . . , n and 8i for i = 3, . . . , n — 1. Now let 8 = 8 n -i 
and let x* be any point of type IX such that d(x*,z n ) < 8. We claim that the solution 
with initial data for r = has the desired properties. Note first that by construction it 
enters each of the balls Bi and is contained in S. Let r n = 0. Let a n -\ = r n + Ar n _i. Let 
r n -i > o~ n -i be a time such that <i(/(r n _i), z n -i) < e n _i and /(r) is in -B^-i for all r in 
the interval [o"j_i, Tj_i]. Continuing in this way, the other Tj can be defined recursively. 

The intuitive meaning of this theorem is that the solution visits small neighbourhoods of 
the points generated by th BKL iteration and, moreover, does so in the order determined 
by the iteration. The analogous theorem holds with type VIII replaced by type IX, the 
proof being almost identical. This result shows that arbitrarily long but finite segments of 
a BKL iteration can be realized by a solution of the Einstein equations. However, it says 
nothing about the full iteration, if it happens to be infinite. This is natural, since for a 
given measure of error e and a given starting point specified with precision e it cannot be 
expected that the qualitative behaviour of the sequence of iterates is determined. 

What do Theorems 5.1-5.3 tell us about Mixmaster dynamics? They tell us that a 
finite sequence of BKL iterates corresponds to a solution of the exact equations. They 
do not tell us the ultimate fate of the solution at early times because the possibility is 
left open that eventually all solutions show the behaviour described above as anomalous. 
Numerical results make this possibility implausible. Recent calculations [4] have been able 
to observe 150 BKL iterates in numerical solutions of the Mixmaster model. Moreover, 
there is no indication that things change at that point and anomalous behaviour seems 
never to have been observed numerically. Thus it seems reasonable to suppose that at least 
a large open set of initial data for the Mixmaster model gives rise to standard oscillatory 
behaviour. On the other hand, it it difficult to see how the mathematical techniques used 
in this paper could lead to a proof of this fact. The key question is, what other techniques 
might do so. 
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